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Abstract 

We consider second-order uniformly elliptic operators subject to Dirichlet 
boundary conditions. Such operators are considered on a bounded domain 
f2 and on the domain <^(J1) resulting from by means of a bi-Lipschitz map 
(f>. We consider the solutions u and u of the corresponding elliptic equations 
with the same right-hand side / G L 2 (Q U (f)(Q,)). Under certain assumptions 
we estimate the difference \\Vu — ^u\\L 2 (nu<t>(Si)) ^ n t erms of certain measure of 
vicinity of <p to the identity map. For domains within a certain class this pro- 
vides estimates in terms of the Lebesgue measure of the symmetric difference 
of (f>(fl) and fi, that is |0(O)AO|. We provide an example which shows that 
the estimates obtained are in a certain sense sharp. 
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1 Introduction 



Let Q be a bounded domain in M. N and let A = {Aij(x)} be a measurable positive- 
definite symmetric matrix-valued function on M. N bounded away from zero and in- 
finity. Let 

i,j=l ■> 

subject to Dirichlet boundary conditions. Let : Q — > be a bi-Lipschitz map 

and let L be the analogous operator on L 2 (<f)(Q)). We fix / G L 2 (Q U <f>(£l)) and 
consider the functions u G Hq(Q) and u G H^(<f)(Q)) denned by 

Lu = f, in Q , Lu = f , in <ft(Q). 
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Our aim in this article is to estimate ||Vtt — Vw||i2(nu0(n)) in terms of a certain 
measure of vicinity of to the identity map. 

Domain perturbation problems are an important branch of the theory of PDEs. 
Within spectral theory one is typically intersted on the stability of the eigenvalues 
or eigenf unctions of differential operators. There are several recent results on this 
type of problems; we refer to the article [7] and references therein for more on recent 
progress on domain perturbation problems in spectral theory. 

The problem we are interested in involves the stability of the solution u of the 
equation Lu = f under perturbation of the domain Q. In the article [9] Savare and 
Schimperna obtain very detailed sharp estimates on the variation in Hq of u for 
operators with Lipschitz continuous coefficients. In [4] estimates are obtained on the 
variation of u in L 2 for operators with measurable coefficients. See also the articles 
[U El [5] where relevant results were obtained. 

In the present article we consider uniformly elliptic operators with measurable co- 
efficients and we prove stability estimates in Hq for the solution u of Lu = f. A 
simple example shows that the estimates are in a certain sense sharp. Our main 
assumptions are, roughly, that Q is perturbed by a global bi-Lipschitz map and 
that / and Vw belong in L q for some q > 2. We note that if / G L q and Q has 
Lipschitz boundary in the sense that there exists a bi-Lipschitz transformation that 
maps Q onto a domain with C 1 boundary, then Vu G L q by a well known result of 
Meyers [8], so our results are applicable. The proof relies on the so-called pull-back 
method. The operator L that acts on <fi(Q) induces naturally an operator on Q, and 
it is that operator which is then compared with H. Hence, in the first section we 
prove a stability estimate for operators acting on the same space under variation of 
the coefficients. 

The method of proof easily generalizes to other kinds of differential operators, such 
as higher-order operators, operators subject to Neumann boundary conditions or 
operators acting on Riemannan manifolds. For the sake of simplicity and brevity 
we restrict our attention to second-order Dirichlet operators on bounded Euclidean 
domains. 



2 A general stability theorem 

In this section we prove an auxiliary result which we believe is of independent interest. 
Let Q be a bounded domain in M. N and for e > let A e = (A l J) be a family of real, 
symmetric, matrix- valued measurable functions on Q satisfying 




(1) 
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For e > we define the self-adjoint operator 




on L 2 (Q), 



subject to Dirichlet boundary conditions on dQ. We now fix / e L 2 (Q) and we 
denote by u e the solution of H e u e = f, e > 0. Using a standard argument we obtain 
from ([I]) that ||w e ||fri(n) < c for all e > hence, up to a subsequence, (u e ) converges 
weakly in Hq(Q). 

We now assume that \\A e — A)||lp(c) = sup ij {||A* jf — Aq \\lp(0)} for some (equiv- 
alently, for all) 1 < p < oo. We then easily deduce that the weak i7 1 -limit of (u t ) 
is precisely u . This implies in particular that the full sequence (u e ) converges to 
no weakly in Hq(Q). It is a natural question to seek conditions under which the 
convergence u € — > uq is strong in Hq(Q). To our knowledge this problem has not 
been studied. In the next theorem we provide four conditions each one of which 
guarantees strong convergence in Hq(Q). The first of these conditions also provides 
an estimate for the rate of convergence u e — > u in Hq(Q), and is the one that will 
be used for the domain perturbation problem in the next section. 

Given a real symmetric matrix A, the matrix A + is defined by means of the spectral 
theorem: if A = A n (e n <g> e n ) is the spectral representation of A, then A + : = 
^2(\ n ) + (e n <S> e n ). 

Theorem 1 Assume that any one of the following four conditions is satisfied: 

(1) Vm G L q (Q) for some q > 2. 

(2) The eig en functions {4> n } of H Q satisfy ||V0 n ||i«(o) < cXI for 
some q > 2, 7 > and all n£N. 

(3) There exists a compact set K C Q such that A e — > Aq in L°°(Q \ K). 

(4) (A -A e ) + -^0 inL°°(Q). 

Then u e — > uq in Hq(Q) as e — > 0. Moreover, in case (1) we have the estimate 

\\Vu e - Vu \\ L 2 (n) < c||Vtto||L«(n)||oe - a \\ L 2 q /( q -2) {n) , e > 0. (2) 



Proof of TheoremUi Part (1). Using the standard notation of repeated indices, we 
have 

tjduo ( due_ _ duo )d _ [ f( )d 

OX % OXj OXj Jq 
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Subtracting yields 

^ _ ft* j*. _ ft* = /• ^ _ ^^p. _ p )dx _ (3) 
n raj axj OTj ax, axj ax, ax,- 

hence, by (JT]) and Holder's inequality we obtain 

||Vu £ - Vw |U 2 (Q) < c||Vu |U«(n)||A 6 - ^o|U%/(8-2)(n)- 

Part (2). Let us denote by V the gradient operator, from Hq{VL) to L 2 {Vt). By 
[31 Lemma 3] we have 

u-u Q = V*A\l\A\/ 2 VV*A\l 2 +lY 1 A^I\A Q -A e )A- l/ \Al / \v*Al /2 +iy 1 Al /2 Vf, 
therefore 

||Vu e - Vu \\ L 2 {n) < \\A: l ' 2 \\ ■ \\A 1 J 2 VV*A l J 2 (A 1 J 2 VV*A 1 J 2 + l)- 1 )] x 

\\A- e l l 2 {A - A e )A- 1/2 (Al /2 VV*A l J 2 + 1)^^ /2 V|| • ||/|| i2(a) 
< c\\(A - A e )A- 1/2 (Al /2 WV*A l /2 + l^A^VW , 



where non-indexed norms are operator norms on L 2 (Q). By [3], Lemma 4] we conclude 
that 

||Vw e - Vu \\ L 2 (n) < c\\A e - A \\ L r (n) , all r > -^(^ + 2 7 - 1), (4) 

q — 2 2 

and the result follows. 

Part (3). Going back to ((3]), applying ([T|) to the left-hand side and decomposing the 
integral in the right-hand side by integrating in K and in Q \ K, we get by Holder 
inequality for any q > 2, 

^||Vu e - Vu ||| 2(n) < C\\A t - A \\ L oo {n \ K) \\Vu e - Vu \\ L 2 (n) + 
C\\A e - A \\ L 2 q /( q -2 ){I q\\Vu \\ Lq{K) \\Vu t - Vu \\ L 2 {n) 

which implies, 

||Vw e - Vn || L 2(n) < C(\\A e - A \\ L oo {n \ K) + \\A e - A \\ L 2 q /( q -2)^\\Vu Q \\ L i(K)) 

But interior estimates for the limit problem imply, see [SI Theorem 2], that there 
exists a q > 2 such that || Vwo||i?(K) < C. Hence, we obtain the stated result. 

Part (4)- Let us start noticing that in general, since u e — > u weakly in H 1 ^) and 
in particular in L 2 (Q), we get 

^du.du, = f fUe ^ f fuQ= f ^^fo^o (5) 
e dxi dxj J n J n J n dxi dxj 
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Hence we have 



J n dxi dxj J n dxi dxj J n dx { dxj 

— > / An '-?r- --^ — dx as e — > 0. 

where we have used the hypothesis and (jSJ). Hence 

f ijdu e du e 2 f ijduoduo 2 

hmsup / (A/ — h tt £ )dx < / (A 7 -— h u )dz. 

e^o J n dxi dxj 

Let us now consider the space H l (Q) with the Hilbert norm 

We have thus proved that limsup |||u e ||| < |||%|||- Moreover, the weak convergence in 
Hq(Q) with respect to the standard norm implies the weak convergence with respect 
to the equivalent norm |||-|||, hence liminf |||u e ||| > |||%|||. We thus conclude that 
|||w e ||| — > |||ito|||) an d therefore u e — > uq strongly in Hq{VL), as required. □ 

Remark. If Q has a Lipschitz boundary, in the sense that there exists a bi-Lipschitz 
transformation that maps Q onto a domain with C l boundary, then condition (2) is 
satisfied. This follows easily from [8j Theorem 1]. 

Remark. It is an interesting problem to prove or disprove that u t — > Uq in Hq(Q) 
without any assumptions other than those stated before Theorem [TJ 

We now present an example that shows that estimate ([2]) is sharp. We fix an angle 
7T < f3 < 2n and denote by VLp the circular sector of radius one and angle (3, 

= {(r,6) : < r < 1 , 0<6 <(3}. 

Let Mo be the solution of the problem, 

-Au = ^sin(f), infy, 

Mo = 0, on dVLp. 

(The factor (4/3 2 — n 2 )/f3 2 is introduced to simplify subsequent calculations.) Simple 
computations give 

u (x) = (r* -r 2 )sin(^) , < r < 1, < 9 < (3 . (7) 

This implies in particular that V«o £ L q (Vtp) if and only if q < 2/3/ (/3 — ir). We now 
fix a > and for < e < 1 we set 



<xax) 



[ a, 


< \x\ 


U 


€ < \x\ 



We note that a e — > 1 in L p {VLp), 1 < p < oo, but not in L 00 ^^). Let us denote by 
u t the solution of the perturbed boundary-value problem 



-div(a e V« e ) = sin(f ), in Q p , 

u e = 0, on dflp. 

From part (2) of Theorem [I] follows that 

2 

||Vw e - VM ||L2(n^) < c\\a e - a \\LP(n p ) = ce? , all p > 
Now, simple computations give that u e (x) = v e (r) sm(ir9//3), where 



(8) 



2(3 



7T 



v e (r) 



( e-l)(^ 2 + ^) + 2e-f f 1 2 
(1 - a)e% + (1 + a)e~% a 



< r < e, 



(1 - a)e 2 + (1 + a)e f> * 



(1 - a)ee + (1 + a)e ? (1 - a)e^ + (1 + a)e 
Let us denote «o(x) = vo(r) sm(n9/ (3) (cf. (0)). We then have 

||Vu e - Vuollia(n-) > c [ (v' e -v' ) 2 rdr 

Jo 



(l-a)(e?-e 2 ) _ s 2 

— r f> - r , e < r < 1. 



7r(l — a) 



,6^(1 +0(1)) 



as e 



0. 



2/3(1 + a) 2 ' 

We conclude that the index 2q/ (<? — 2) in (j2J) cannot be replaced by any smaller index. 



3 Stability estimates in Hq under domain pertur- 
bation 

Let n be a bounded domain in R N . We fix M > 0. Let : Q -»■ 0(0) =: O be a 
bi-Lipschitz transformation such that 

IWIU~(n) < M, IK^)- 1 !!^^ < M . (9) 
We set E = {x E Q : <f>(x) ^ x} and we note that 

\\F o <f) - F\\ L 2 {n) < ||F|| L9(n) |£|^ , IKD^-Jll^n) <c|E|^. (10) 
Let A = {Ajj(x)} be measurable and symmetric on Q U Q satisfying 

^\£\ 2 <Y^Mx)tej<M\£\ 2 , xeQun, £eR N . 
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Let L (resp. L) denote the operator — ^2ijd Xj {Aijd Xi } on L 2 (Q) (resp. L 2 (Q)) 
subject to Dirichlet boundary conditions. We fix / G L 2 (Q U f2) and we denote by u 
(resp. u) the solution of Lu = f (resp. Lu = f). We extend Vm and V£t to be zero 
outside their respective domains. We then have the following 

Theorem 2 Assume that there exists q > 2 such that H/Hj^mun) < M, ||Vw||i«(n) < 
M and ||Vn|| L ,^^ < M. T/ien there exists a constant c depending only on M such 
that 

\\Wu-Vu\\ L2(nu ^<c\E\ q -W. (11) 



Proof. We set g(x) = | det D(j){x)\, igfl. We note that given v G Hq(Q) and letting 
u = v o we have 

/ |f | 2 <iy = / |u| 2 g ofo; 
</n Jn 

and 

where a = (aij)i,i=i,-,JV ^ s defined on f2 by 

^ rj i ( — 1) O I ( — 1) 

£ (A rs -| 1*—) o = ( W)" 1 ^ o <f>){D<f>)-% . (12) 



r,s=l 



This leads to the notion of the pull-back: the pull-back of L to fl is the self-adjoint 
operator H on L 2 (Q,gdx) associated with the sesquilinear form Q with Dom(Q) = 
H£(tt) and 

~ f J^, duidu 2 _ rr] 

Q{u 1 ,u 2 )= yaij — —gdx, u 1 ,u 2 eH (tt). 

Jn i j=1 ax i ax j 

So formally if?7 = —g~ ls 22ii j{ga>ijU Xi ) Xj . Equivalently, if can be defined as H = 

CfiLC^ 1 , where : L 2 (fl) — > L 2 (Q,gdx) denotes composition by (a unitary 
operator). 

It then follows that H = g _1 if, where H is the self-adjoint operator on L 2 (Q) 
associated to the form Q defined above on Hq(Q). Now, let u G L 2 (Q) be defined by 
Hit = f. Using (1) of Theorem [1] and also ffTUl) we obtain 

II V« — Vnll r2/ni < ell a — All 2q 
ii ii^ W - II M L^(n) 

< c\E\~*r. 
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Morover we have 

H(u o <f>) = / o 
Hu = gf. 



(13) 



Hence, by (ITO 



|V(wo0)- Vu|| L 3 ( n) < c||/o0-^/|| i2(n) < c|£7|V. (14) 



Finally by ( llOp we also have 



q-2 



\V(uo(j))-Vu\\ L 2 {n) <c\E\ — . (15) 



3-2 



Combining (|T5|) . ( TClj) and ([15]) we obtain that ||Vu - Vw|| L 2 (f7) < c\E[W . We also 
have 

< c|0 -1 (fi\fi)| < c|£|, 

q-2 

hence ||Vw — Vw|| L 2(^) < c l^| 2<? by Holder inequality. This concludes the proof. □ 

Using Theorem |2] we can now prove stability estimates in for localized perturba- 
tions in terms of the Lebesgue measure. We consider the following class of domains 
(see also [6]): 



Definition 3 Let V be a bounded open cylinder, i.e., there exists a rotation R such 
that R(V) = Wx]a, b[, where W is a bounded convex open set in IR^ -1 . Let M, p > 0. 
We say that a bounded open set Q C W N belongs to (V, R) ifVt is of class C 0,1 (i.e., 
Q is locally a subgraph of C 0,1 functions) and there exists a function h G C°' l (W) 
such that a + (b — a)/10 < h < b, Lip(/i) < M, and 

R(Q n V) = {(x, x N ) : x eW , a < x N < h(x)}. (16) 



In the following theorem, again, V« and Vw are extended to be zero outside O and 
fl respectively. 

Theorem 4 Let f2, Cl be domains of class C^(V, R) with Q\V = Cl\V . Let f G 
L 2 (Q U Q) and let u G Hq(Q) and u G Hq(Q) solve Lu = f and Lu = f respectively. 
Assume that there exists q > 2 such that ||/||i«(nun) — M, || Vu||l«(J2) < M and 
||Vw|| w ^ < M. Then there exists a constant c depending only on M such that 

l|V« - Vu\\ L2(auCi) < c\QAQ\^ . (17) 
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Proof. It has been proved in [I], following earlier work of Burenkov and Lamberti 
[6], that if f2 and Cl are in (V, R) then there exists a bi-Lipschitz diffeomorphism 
(j) : Q, — >■ Cl whose Lipschitz constants are estimated in terms of M only and such that 
\E\ < c\ClAQ\, where c also depends only on M. Hence (j!7p follows from Theorem 

m □ 

We now present an example that shows that estimate ([2]) is sharp. It is a variation 
of the example in Section [2j We fix ir < f3 < 2ir and < e < 1 and denote by Vt e ^ 
the planar domain 

tt e , p = {(r,0) : e < r < 1 , < < /?}. 
We denote by u e be the solution of the problem 

-A« e = ^sin(f), infi^, 

(18) 

u e = 0, on dVt e ^. 

The function uq has been computed in ([7]). For e > a direct computation gives 
that u t (x) = v e (r) sin(7T0//3), where 

v e {r) — jr. t*tp H -r, r--r 0— r , e < r < 1. 

We then easily obtain the asymptotic formula 

||Vu e - Vuollrw = ^(1 + o(l)), (19) 

for some A > 0. 

Let us now see what our theorem gives. We note that the domains are not of some 
class Cm~(V,R) uniformly in e, so we cannot directly apply Theorem HI Instead, let 
s e :(0,l)^(e,l), 



sJr) 



| + e, 0<r<2e, 



r, 2e < r < 1. 
We define the bi-Lipschitz map e : Qq — > f2 e , given in polar coordinates by 

e (r, 0) = (s e (r), 0) , < r < 1 , < < /3. 

Using the explicit computation of Wo an d w e we easily see that the assumptions of 
Theorem [2] are satisfied for any (fixed) q < 2/3/(/5 — tt). Hence Theorem [2] gives 

||V« e - VM ||L 2 (n un E ) < Cg|-E e |^" = c ? eT = c 5 e^" 5 , 

for any 5 > 0. Because of fTT9l . this shows that the exponent of \E\ in fTTTl) cannot 
be replaced by a smaller one. 
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